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We study the asymmetry between the production of baryons and that of antibaryons in relativistic heavy ion 
collisions in quark combination model. Based on basic ideas of quark combination, we study the general features 
of the bulk production of baryons and mesons in hot and dense system and obtain yields of baryons, antibaryons 
and mesons as functions of the numbers of constituent quarks and antiquarks just before hadronization. In 
particular, we study effects of quark-antiquark asymmetry on the ratio of (anti-)baryons to mesons and that of 
antibaryons to baryons. The universal correlation between pip and K~ /K + yield ratios observed in heavy ion 
collisions at energies ranging from SPS to LHC is explained. 
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I. INTRODUCTION 

Heavy-ion collisions at relativistic energies produce a 
long lifetime bulk quark matter with extremely high tem- 
perature and energy density yj. The system has a nonzero 
baryon number density due to the ever present nuclear stop- 
ping J2l- This leads to the asymmetry between the produc- 
tion of hadrons and that of antihadrons, measured mainly by 
the ratios of yields of antihadrons to hadrons [3]. Production 
of baryons and antibaryons, compared with that of mesons, 
is strongly influenced by the baryon number density of the 
system due to the baryon number conservation[3]. 

Ratios of yields of antiprotons to protons and that of K 
to K + are two of the most significant observables measuring 
the hadron-antihadron asymmetry in heavy ion collisions. 
Final state proton ratio p/p carries the information of the 
production of baryons and antibaryons. Kaon ratio K /K + 
roughly cancels the effect of strangeness production and re- 
flects the asymmetry between charged mesons and their an- 
tiparticles. Experimental data of this pair of baryon and 
meson ratios are available from the Relativistic Heavy Ion 
Collider (RHIC) from the relatively low energy col- 

lisions at the Super Proton Synchrotron (SPS) and Alternat- 
ing Gradient Synchrotron (AGS) Il8l-ll0ll and more recently 
from the very high energy reactions at the Large Hadron 
Collider (LHC) ll ill . The data seem to suggest a universal 
correlation between the p/ p ratio and K/K + ratio at these 
collisional energies [5]. Popular event generators A Mul- 
tiphase Transport Model (AMPT) default version lfT2ll and 
Ultra-Relativistic Quantum Molecular Dynamics (UrQMD) 
lfl3ll . both dealing hadronization with string fragmentation, 
cannot explain this phenomenon j^. We note that the re- 
sults of these two generators are similar to that of PYTHIA 
for pp collisions. This indicates that hadronic rescattering 
effects in heavy ion collisions have little relevance to the 
observed correlation which should be closely related to the 
hadronization of bulk quark matter. 

In this paper, we use quark combination mechanism 
(QCM) to explain the observed correlation between p/p 
and K /K + ratio. Quark combination attributes the ratios 
of antihadrons to hadrons to the asymmetry between the 
number of constituent quarks and that of antiquarks before 
hadronization. K /K + ratio is mainly determined by N q /N q 
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because of strangeness neutrality N s = TYj while p/p ratio 

depends on the high power of N q /N q , e.g. p/p ~ [Nq/Nq) 
in a naive quark counting. Here, q and q denote light quark 
and antiquark, respectively. In this view, the final state p/p 
ratio as the function of K /K + ratio has an intuitive physical 
meaning — it reveals the law of the production asymmetry 
between baryons and antibaryons with respect to the baryon 
number density of the system characterized by the quark- 
antiquark asymmetry. Therefore, we will emphasize on the 
bulk production of baryons and antibaryons in QCM. 

The paper is organized as follows. In Sec II, we will start 
from the basic ideas of the combination mechanism to carry 
out a general analysis of the yield of baryons and that of 
mesons. We study the dependence of yields of directly pro- 
duced baryons and mesons at hadronization on the number 
of constituent quarks and antiquarks. In particular, we are 
interested in the influence of quark-antiquark asymmetry on 
the yield ratio of (anti-)baryons to mesons and that of an- 
tibaryons to baryons. We will make the study as indepen- 
dent of the particular models as possible but present the as- 
sumptions and/or inputs explicitly whenever necessary. In 
Sec III, we make the comparison with the experimental data 
of the p/p ratio with respect to K /K + ratio in heavy ion 
collisions at energies ranging from SPS to LHC. Sec IV 
summaries our work. 



H. GENERAL ANALYSIS OF HADRONIC YIELDS IN 
QCM 

Considering a bulk quark system consisting of N q con- 
stituent quarks and N q antiquarks, after hadronization the 
system changes the basic degrees of the freedom to become 
the hadronic matter and produces in average 

B(N q ,N q ) baiyons, 

B(N q ,N q ) antibaryons, 

M(N q ,N q ) mesons 

which are functions of N q and Nq. Here we temporarily 
leave the hierarchy properties of quarks such as flavor and 
mass aside to focus on the global production of baryons and 
mesons. M and B(B) refer to the number of all mesons and 
(anti-)baryons, respectively. 

An important property of yield functions of hadrons 
comes from the charge conjugation transformation of quark 
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to antiquark. Under this transformation, baryons become 
antibaryons, antibaryons become baryons and mesons still 
keep mesons. Hadron yields change as B ^ B and M ±=t M 
under q ^ q exchange, and thus we have 



B(N q ,Nq) = B(Nq,N q ), 
M(N q ,N q ) - M(N q ,N q ). 



(1) 



Quark combination mechanism describes the production 
of hadrons by the combination of constituent quarks and an- 
tiquarks, that is, a quark and an antiquark merge into a me- 
son and three quarks into a baryon. After system hadroniza- 
tion, there are no free quarks and antiquarks left and all of 
them are combined into hadrons. We use the following con- 
stituent quark number conservation to denote the finish of 
hadronization process, 



M(N q ,Nq) + 3B(Nq,Nq) = Nq, 

M(N q ,Nq) + 3B(N q ,Nq) = N q . 



(2) 



This unitarity is often adopted in studying hadronic yields 
in literatures [ 14-16]. The related issues of entropy and en- 
ergy conservation have been properly addressed in [fnl HH 
by considering volume expansion and temperature decrease 
during system hadronization. Here we neglect the produc- 
tion of multiquark states and exotic hadrons. Since we are 
interested in hadronic yields which are low pj dominated, 
we also neglect the hadronization of jet quarks and the inter- 
action of jet quarks and thermal quarks in the matter because 
of their negligible contributions to hadronic yield. 

In addition, as the system entirely consists of quarks or of 
antiquarks, quark combination, if still works, should give 



\B = 0, 

N: 



N 

B = !f,M 







if Nq 







B=-f,B = 0, M = if N q = 



(3) 



These are boundary conditions of yield function of baryons 
and mesons in quark combination. 

Obviously, the yield of baryons/antibaryons under bound- 
ary conditions Eq.Q is the linear function of that of 
quarks/antiquarks. In general case of both mesons and 
baryons being produced, supposing that we synchronously 
increase (or decrease) the number of quarks and that of 
antiquarks to vary the size of quark system at hadroniza- 
tion without changing the intensive properties of the system 
such as temperature and baryon number density, intuitively, 
hadronic yields should vary with the corresponding magni- 
tude, i.e. linearity holds also. This linearity is supported by 
the RHIC data which show that the multiplicity of charged 
particles in central and semi-central nucleus-nucleus colli- 
sions is of nearly linear proportion to the number of partici- 
pant nucleons jV garf fl9tl while the pln~ yield ratio, even at 
LHC energies ifllh - is almost independent of N part fK. l20ll . 
Calculations of the quark combination model developed by 
Shandong Group (SDQCM) [21] also support such a prop- 
erty. Letting A is a not very small auxiliary parameter, we 
thus have 



M(AN q , ANq) = AM(N q , Nq), 
B(AN q , ANq) = AB(N q , Nq), 

B(AN q , AN q ) = AB(N q ,N q ), 



(4) 



in quark combination hadronization. Note that here we 
focus on the bulk matter with abundant particles and do 



not consider the situation of small system with few quarks 
where extra restraints in production of hadrons, in particular 
of baryons and anti-baryons, maybe come out. 

These properties Eqs. ([T]|4j pose important constraints on 
the behaviors of yields of baryons and mesons. First, we 
rewrite the arguments of yield function, i.e., N q and Nq, in 
terms of x and z by 



N Q + N s = x, 



N q -N q 
N q +N q 



= z. 



(5) 



Variable x is the total number of quarks and antiquarks, 
which characterizes the bulk property of the system related 
to the system size or energy. Variable z depicts the asymme- 
try between quarks and antiquarks in the system (\z\ < 1), 
which is a measure of baryon number density of the system. 
Under q ^ q permutation x ^ x and z —z and yields of 
baryons and mesons can be written as the function of x and 
z 

B(N q , Nq) = B(x,z), (6) 
B(Nq,Nq) = B(x,z) = B(x,-z), (7) 
M(N q ,Nq) = M(x,z) = M(x,-z), (8) 

where we have used the symmetry properties Eq. ([1}. Cor- 
respondingly the linearity property Eq. (0) can be expressed 
with variables x, z as 



M(Ax,z) = AM(x,z), 
B(Ax,z) = AB(x,z), 
B(Ax,z) = AB(x,z). 



(9) 



Since the B(x, z) is the linear function of variable x, we im- 
mediately have 

B{x,z) = xS(z), (10) 
and correspondingly for antibaryons and mesons, 

B(x,z) = B(x,-z) = xS(-z), (11) 
M(x,z) = M(x,-z) = |jc(l -3S(z)-3S(-z)). (12) 

Properties of unitarity and boundary conditions at 
hadronization exert constraints on S(z) as follows 

!B(l) = \, £(~D = 0, S(z)~!B(-z) = z/3. (13) 

Using Eqs. ( ITOb and ( fT3l we obtain the ratio of yield of 
anti-baryons to that of baryons. 



B(x,z) S(-z) . 
R(z) = — r = = 1 



B(x,z) S(z) 
The ratio R satisfies 

*(0)=1, 
K(1) = 0, 
R(-l) = 00, 

R(z)R(-z) = 1. 



3S(z) 



(14) 



(15) 



We can see that the behavior of R(z) is greatly constrained by 
properties Eqs.(Q]|4]i which come from the global properties 
and/or symmetries. 
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The complete determination of R(z) needs further con- 
sideration of hadronization dynamics. In QCM, baryon is 
formed by the combination of three (anti-)quarks in neigh- 
borhood in the phase space. The (anti-)baryon yield is de- 
termined by two ingredients: (1) the probability of finding 
three neighboring (anti-)quarks in phase space P(qqq) mq . (2) 
the probability of these three (anti)quarks to form a (anti- 
)baryon, P( qqq ^,B), qqq ^B- The contribution of (1) to R(z), 
Pqqq/Pqqq, follows the form of (N q /N q Y for the system 
consists of free quarks and antiquarks. The combination 
probability P( qqq -,B), qqq ^B m (2) is determined by the local 
hadronization dynamics. If the fast hadronization approxi- 
mation is assumed as did by most of combination models, 
the P qqq ->B is evaluated by the overlap between the hadronic 
wave function and quark ones, which is undisturbed by sur- 
rounding environment. Then it contributes nothing to R(z) 
and we can obtain R(z) = (N q /N q ) a . On the other hand, if 
we consider the possible effects of the surrounding environ- 
ment (i.e. quarks and antiquarks) to qqq — > B process, we 
may introduce a further contribution to R(z) because the sur- 
rounding quarks and antiquark are asymmetrical. Supposed 
that there is a q in the neighbor of qqq cluster, this q can 
have the chance to capture one q to form a meson and left 
two quarks combining with other quarks/antiquarks. On the 
other hand, if the nearest neighbor of qqq cluster is a quark, 
the baryon formation should be then dominated. Therefore, 
the ingredient (2) can also contribute to R(z) via the form 
of Nq/Nq. We note that such a consideration is essentially 
a kind of dynamical competition between the production of 
baryon and that of meson. Finally, we obtain 
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/ 1 — 



(16) 



by noticing the identity N q /N q — (1 -z)/(l +z). It is the sim- 
plest functional form with one parameter a which satisfies 
Eq. ( TT5b . This form can be tested by an intuitive combina- 
tion process in the following way. First, we sort all quarks 
and antiquarks in the system into a one dimensional quark 
sequence in which the position of quarks can denote the de- 
gree of proximity to each other. Second, we combine those 
qq and qqq in neighbors into mesons and baryons, respec- 
tively, until all quarks and antiquarks are exhausted. 

Parameter a bears important hadronization dynamics 
which is reflected by 

3^b/m(0) 

where R B /m(0) = B(x,0)/M(x,0) is the ratio of the yield 
of baryons to that of mesons at vanishing quark antiquark 
asymmetry z = 0. This ratio denotes the competition be- 
tween the formation of baryon and that of meson when 
quarks are combined into a hadron. The study of hadron 
yields in Ref. ir22ll shows that Rb/m(0) of value about 1/12 
can well describe the data of p/V~ yield ratio in Pb+Pb col- 
lisions at = 2.76 TeV after including effects of reso- 
nance decay. Then we obtain a * 5. 

The functional forms of B(x, z) and B(x, z) then read as 

(1 + zY 



B(x,z)=x!B(z) = 



B{x,z) = x!B(-z) 



xz 

~'(l+zY 

xz (1 



(1 -zY' 
■z) a 



3 (1 +z) a -(1 -zY 



(17) 



(18) 



Figure Q] shows the results of Rb/m(z) and R b /m(z) = 
Rbim(-z). Symbols are numerical results of SDQCM yJJ. 



B/M SDQCM 
□ B/M SDQCM 




FIG. 1: (Color online) The ratio of yield of (anti-)baryons to that 
of mesons as the function of z. 



III. COMPARISON WITH EXPERIMENTAL DATA 

In this section, we test our baryon formula by the avail- 
able experimental data of pip at RHIC and SPS. Since quark 
asymmetry z is unobservable we should firstly relate it to 
an observation, Nk-/Nk+ yield ratio. Now we consider that 
(anti-)quarks have three different flavors, i.e. u, d and s. 
Their numbers are denoted by N u ^,N d j andA^j. Under 
isospin symmetry we can write z as 



N q -N q 
N q + N q 



N u -N Tl 
N u +N a +N s 



1 _ !k 

N„ 

1 + ^ + ^' 

N„ N„ 



(19) 



In the quark combination mechanism, ratios Nu/N u and 
N s /N u can be related to the corresponding ratio of 
different hadron species that contain these constituent 
quarks. For example, we have Nk-(us)/Nk+( U s) k N b /N u 
and N,f,( S s)/NK+(us) x N s /N u in general because of the 
strangeness neutrality N s = N s - Adopting our terminology 
of the description of quark combination in Ref. ll2.2ll for the 
proportion coefficients, these two ratios read as 

Nk_ = C^Nn^ = 

N K * Ck+Nu N k + C k+ N u ' 

Here, Ck* and are combination factors denoting the 
branch ratio of a given constituent quark composition (e.g., 
us) to form a specific hadron species (e.g., K + ). Ck- = Ck* 
due to both hadrons are in the same J p = (T pseudoscalar 
nonet. C$/Ck+ = Rvp where Ryp represents the ratio of the 
J p — 1~ vector mesons to the J p = 0~ pseudoscalar mesons 
of the same flavor composition. 

The influence of resonance decay on Nk-/Nk* ratio is 
small because of almost same decay contributions to K ± , 
and we still have 



(final) 



N (J'inal) 



(21) 



For N^/Nic* ratio, decay of resonances contributes signifi- 
cantly the denominator, and thus we have 



N 



(final) 



R 



VP 



(final) 



N< 



N K * +N k , 1 + R V /p N u 



(22) 
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We further rewrite l.h.s of the equation as 



,v 



by 



noticing that the experimental data of jv^ m£ ")/jv^. is 
an almost collision-energy independent value, which is 
about 0.14 |[23tl . Then we can relate the argument z to 

N (fmal) /N (fmal) (denote( j by K ~ j K + fo r short ) as following 



z = 



1 - £- 

1 K+ 



1 + 



(l +0 .14^)§ 



(23) 



The variance of /?y/p exerts little influence on the transition 
between z and K /K + . Here we take Ry/p = 0.45 based on 
the study of K* production in relativistic heavy ion collisions 
inRef.Q. 

By definition, R(z) is the ratio of the yield of directly- 
produced anti-baryons to that of baryons in the system at 
hadronization. Since all (anti-)baryons will decay into (anti- 
)nucleons at last, R(z) equals the ratio of final state anti- 
nucleons to that of nucleons and finally we obtain 



R( Z ) = N ^ a 'l N f p inal 



(24) 



under the isospin symmetry, which is denoted by p/ p for 
short. Now, with Eq. (1231 we can use the data of K /K + to 
predict the final state p/p ratio to test our formula of baryon 
production and underlying assumptions. 

Figure |2] shows the p/p ratio with respect to K/K + ra- 
tio in central nucleus-nucleus collisions at energies rang- 
ing from ^JTNN - 5 GeV - 2.76 TeV. The data points 
of RHIC energies JH-01 are obtained in different rapidity 
slices, while the SPS and AGS points I18l- 4l Qfl are obtained 
at midrapidity. The data of LHC in central Pb+Pb collisions 
at -^Sffff = 2.76 TeV 111 ill are also presented. Comparing the 
low and high energy results, we observe a coincidence of the 
forward rapidity results of Au+Au 62.4 GeV and those at 
SPS energies, and also a coincidence of the midrapidity re- 
sults of Au+Au 62.4 GeV and those of Au+Au 200 GeV at 
forward rapidities. This indicates the universality of hadron 
production in heavy ion collisions at these energies. Results 
of event generators AMPT 1.11 with string fragmentation 
and UrQMD 2.3 in Fig. [2] are taken from Ref. |H. They 
give an qualitative description for p/p variance with respect 
to K /K + ratio but a bad explanation for the universality of 
the correlation. We note that results of these two generators 
are similar to that of PYTHIA for pp collisions, which in- 
dicates that this systematic correlation should be relevant to 
the hadronization of bulk quark matter instead of hadronic 
rescattering effects in heavy ion collisions. 

The solid line in Fig. [2] is our result of p/p « R(z) by 
Eq. (TTBT ) in terms of z by Eq. 1231 . We observe a good 
agreement with experimental data. The dashed line is the 
result of SDQCM which shows the influence of full treat- 
ment of resonance decay. In our explanation, the exponent 
a = 5 which determines the slope of B/B ratio is the inverse 
of the products of three factors 3 x 2 x S(0) by using the 
derivation of R(z) at z = 0, ^f!r |__ = - 3 g (0) = -2a. The 
3 comes from the number of constituent quarks contained 
in baryon and the 2 from the symmetry property of baryon 
and antibaryon production. The S(0) is the proportion of 
the baryon yield in the hadronization of the bulk quark mat- 
ter with zero baryon number density. They together predict 
the production properties of baryons and antibaryons at non- 
zero Z- 
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FIG. 2: (Color online) p/p yield ratio with respect to K~ /K + ratio 
in central nucleus-nucleus collisions. The symbols of RHIC en- 
ergies are BRAHMS and STAR data from Refs. [3J47J] and lower 
energy data from Refs. I8l 4l0ll . The data of LHC in central Pb+Pb 
collisions at ysiviv = 2.76 TeV are from Ref.[ll]. Results of 
AMPT 1.1 1 and UrQMD 2.3 are taken from Ref. H]. 



We also show the argument of naive quark combination, 
i.e., p/p - (K~/K + ) 3 by simple quark number counting 
K /K + ~ Nu/Nu and p/p ~ (N ri /N u )\ It is obviously devi- 
ated from (located above) the experimental data. This is be- 
cause such an inclusive point of view only considers the free 
combination of (anti-)quarks at the moment of (anti-)baryon 
production without the proper treatment of final fates of 
all other (anti-)quarks that are failing or are not participant 
in baryon formation. The exclusive quark combination, as 
shown by our result, covers this unitarity of hadronization 
to yield the increased exponent a = 5 for p/p = (K~ /K + ) a 
parametrization, which is in good agreement with the data. 
Description of baryon production in string fragmentation is 
significantly different from that in quark combination. It 
gives a quite large exponent a > 8 and resulting curves lie 
below the data in the intermediate K /K + region. 



IV. SUMMARY 

We have studied the correlation between p/p and K /K + 
yield ratios at RHIC and SPS energies in the framework 
of the quark combination mechanism. These two ratios 
can reflect the property of baryon production in relativistic 
heavy ion collisions. Starting from the basic ideas of quark 
combination, we discussed the general properties of bulk 
hadron production in quark combination and obtained yields 
of baryons, antibaryons and mesons as functions of the num- 
bers of constituent quarks and antiquarks before hadroniza- 
tion. It is found that the observed correlation between p/p 
and K~/K + ratios are well explained by our formula. It is 
important to understand the low pj hadron production in rel- 
ativistic heavy ion collisions. We expect the Beam Energy 
Scan program of STAR Collaboration can complete the data 
set of p/p and K /K + ratio in the intermediate z region to 
make a further test of this universality of hadron production 
in relativistic heavy ion collisions. 

We note that these two ratios have been studied in statis- 
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tical models j25il26ll and can be explained well also. This 
kind of models use the chemical potential to explain ratios 
of antihadron to hadron and thereby gain insight into the 
equilibrium information of the hot and dense matter pro- 
duced in relativistic heavy ion collisions. Our explanation 
of the correlation between p/p and K~ /K + ratio is based on 
a microscopic mechanism, which is different from those of 
statistical models. 
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